A note on the metastability of the Ising model: the alternate updating case 
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Abstract 

We study the metastable behavior of the two-dimensional Ising model in the case of an 
alternate updating rule: parallel updating of spins on the even (odd) sublattice are permitted 
at even (odd) times. We show that although the dynamics is different from the Glauber serial 
case the typical exit path from the metastable phase remains the same. 
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There have been many studies of the problem of metastabihty in the Ising model, i.e. the time 
evolution of a system initially in the all minus phase in a small positive external magnetic field. 
The case of Glauber- Metropolis serial dynamics has been discussed in [1], where it was proved that 
the exit from the metastable phase is achieved via the nucleation of a critical droplet and that the 
life time of the metastable phase depends only on the energy of such a critical droplet. The proof, 
which uses the pathwise approach introduced in [2], is based on the analysis of the tendency of 
droplets of pluses in a background of minuses to shrink or to grow. In [3] it was shown that the 
behavior of such droplets does not change if the spin flip rates are modified. 

In this note we consider the case of parallel dynamics in which spins on the even (odd) sublattice 
are simultaneously updated. More precisely: we consider the two-dimensional Ising model on a 
finite rectangle A C with even side length and periodic boundary conditions; for any x e A, 
(t{x) = ±1 denotes the spin variable at that site and to any configuration cr G O := {—1, +1}^ we 
associate the energy 

<x,y> xgA 

where J ^ h > 0. The equilibrium is described by the Gibbs measure z/((t) := e^p{—(5H{a)}/Z 
with P, the inverse of the temperature, a positive real number and Z := J^^^q exp{—PH{rj)} the 
partition function. 

We define the following updating rule: we partition A into its even and odd sublattices, A = 
A^ U A", with A^(°) := {x = (xi,X2) G A : + X2 is even (odd)} and we denote by S^*^") := 
{{xi, ...,Xm} '■ 1 < m < |A|/2, Xi G A^*^°^ for all i = l,...,m} the set of all possible collections of 
sites in A'^^"). Given / G we denote by the configuration obtained by flipping in a the |/| 
spins associated to the sites in /. Let //^(°) be a probability measure on S^^") such that //^(°) (/) > 
for any / G 3'^^"^ and t = 0,1, . . . an integer variable. Finally, at any even (odd) time t a set 
J g ■£e{o) (.]josen at random with probability and the |/| spins associated to sites in / are 

flipped with probability exp{— /3[i7(cr-'^) — H{a)]'^}, where a is the configuration at time t. In other 
words we can say that the evolution of model (1) is described by a Markov chain at with transition 
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probability P{(Jt-\-i = r]\at = = Pai^: fl)i with a = e, o depending on t even or odd, where for any 



Pa(^,r?) := < 



exp{-/3[ii'(^^) - H{£,)]+} if 3/ G such that r/ = 

(2) 

otherwise 



and Pa(C)0 := 1 — ??)• We remark that this dynamics satisfies the detailed balance 

condition with respect to the Gibbs measure, hence this is its unique invariant measure. 

Our dynamics lives in the same configuration space with the same energy landscape as the 
standard single spin flip Metropolis algorithm, but many more connections between different con- 
figurations have been opened [4]. In other words we allow the system to perform jumps forbidden 
in the standard serial case. To be more precise the notions of "communicating states" (allowed 
jumps) and, hence, of connected sets have been changed: we say that cr, 77 G are e(o)-nearest 
neighbors (e(o)-nn) iff 31 G E^(°) such that 77 = cr^; two configurations are nearest neighbors if they 
are e-nn or o-nn, in other words any two configurations when differing only on the evens or odds 
sides are neighbors. A set ^ C O is connected iff for any a,r] e A there exists a sequence (Jq, cr„ 
of configuration such that: ctq = cr, cr„ = 77 and for any i = 0, n — 1, depending on i even or odd, 
o"i and (Tj_|_i are e(o)-nn. 

The fact that new jumps have been allowed could, a priori, modify the metastable character of 
the dynamics and eventually destroy it. In the most extreme case single spin flip paths containing 
steps against the energy drift can be replaced by direct energy favored jumps. Indeed, let us 
consider a E Q and / = {xj : 1 < i < n} E , it is easy to show that 

H{a')-H{a) = J2m^''')-H{^)] ; (3) 
1=1 

one just has to remark that for any i,j = 1, ...,n and i ^ j the two sites Xi and Xj are not nearest 
neighbors. In other words (3) means that from the point of view of the energy a parallel event can 
be substituted by a suitable sequence of single spin flip events. Now, one can easily choose a and 
I = {xi: 1 < i < n} such that H{a^) < H{a) and i?((7^') - H{a) > for some i = 1, . . . ,n. 
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Nevertheless, by following the scheme of [5] and by making a repeated use of (3) we will show 
that the metastable behavior of the system remains unchanged: the first problem to solve is the full 
characterization of the local minima of the hamiltonian. A local minimum is a configuration o" G Q 
such that for any nearest neighbor configuration 77 one has H{a) < H{r]). Although the hamiltonian 
is the standard Ising hamiltonian, the notion of neighboring configurations is dramatically changed, 
so it is not obvious a priori that the local minima are the rectangles of pluses as in the single spin 
flip case. But this is the case, indeed we can prove that a configuration cr is a local minimum for 
the alternate dynamics if and only if it is a local minimum in the single spin flip case. It is easy 
to see that the sufficient condition holds because a pair of nearest neighboring configurations for 
the single spin flip dynamics is a pair of nearest neighboring configuration even in the alternate 
case. The necessity condition: let o" be a local minimum for the single spin flip dynamics and let r] 
a e(o)-nn of a. There exists I = {xj : 1 < i < n} G such that rj = a^; the proof follows from 
(3) and from the fact that any cr^* is a nearest neighbor of a in the single spin flip case. 

The second step in the understanding of the metastable behavior of the model is the description 
of the tendency to shrink or to grow of the local minima. A preliminary definition: given ry G 
if at is the process started at some cto G O we denote by t^^ := {t > : at = r]} the first hitting 
time on rj. Now we state the following lemma on the criticality of rectangles: a rectangle R^^rn is a 
configuration with a rectangular droplet of pluses, with side lengths i and m, plunged in the sea 
of minuses. 

Lemma 1 Let I* := [^] + 1; consider a rectangle Re^m € with I < m, and the Markov 
chain at started at a^ = Re,m ■ Given e > we have: if i < i* then Re,m subcritical, that is 
P(r_i < T+i) 1, and P(exp{/3(£ - l)h - /3e} < t_i < exp{/3(^ - l)h + ^e}) 1. Ifl> t 
then Re^m is supercritical, that is -P(t+i < r_i) ^-^^^ 1^ and P(exp{^(2J — h) — Pe} < r+i < 
exp{P{2J - h) + Pe}) 1. 

To prove the Lemma we need few more definitions: a downhill path is a sequence of configuration 
(To, <Ji , an such that depending on i even or odd cxj and aj+i are e(o)-nn and for any i = 0, n— 1 
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H{(Ji) > H{ai^i); an uphill path is the reverse of a downhill one. The basin of attraction of a local 
minimum a is the set B{a) := {rj E : all the downhill paths starting from rj end in a}, and for 
any ^ C its boundary dA is the collection of configurations r] £ \ A such that there exists 
a E A nearest neighbor of 77. 

Now, we just have to find the minimum of the energy on dB{Rg^^) and use the results in 
Propositions 3.4 and 3.7 of [5]. Starting from Re^m one should consider all the possible uphill paths 
reaching dB^R^^jn), but (3) implies that all the mechanisms involving more than a single spin flip 
are energetically less favorable; hence the minimum of the energy on dB{R^^m) can be found using 
the same arguments as in the single spin flip case. From this remark the proof of the Lemma 
follows. 

Now we state the theorem which describes the exit from the metastable phase: let us denote by 
V the configuration with all the spins minus excepted those in a x (£* — 1) rectangle and in a unit 
protuberance attached to one of its longest sides. We set T := H{'P)—H{—1) = AJi*—h{£*'^—£*+l), 
we consider the process at starting from —1 and we define f_i := sup{i < r+i : at = — i}, the last 
time the system visits —1 before reaching +1, and f-p := infjt > r_i : at = V}, the first time the 
system reaches V after having "definitively" left —1. Finally, we can state the following 

Theorem 1 Let us suppose ctq = — L let e > 0, then P{f-p < r+i) ^ll^ 1 and P(exp{/3r — ^e} < 
T+i<exp{pr + Pe}'^^ 1. 

We sketch the proof of the Theorem. We use the same notation as in Section 3.1 of [6]. First of all 
wc consider the map : 17 — )■ morally associating to each configuration cr G the "largest" local 
minimum S{a) to which a is connected by a downhill path, where downhill must be understood in 
the sense of one spin nearest neighboring configurations (configurations differing for the value of a 
single spin). See Section 3.1 of [6] for the rigorous definition of S. Now, we define the set „4 C 
as the collection of configurations cr G O such that all the plus spins of S{a) are those associated 
to the sites inside a collection of rectangles on the lattice such that however two of such rectangles 
are chosen, their mutual distance is greater or equal to Vb (pairwaise non-interacting rectangles). 
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Following [5], we prove that H{V) is the minimum of the energy on the boundary of A. Let us 
partition dA into two parts dA = diALid>2A: for each 77 G dA, r) G diA iff there exists x G A such 
that G A, while r] G d>2A iff G r2\^ for all x G A and there exists / G with |/| > 2 such 
that G A. We first prove that for all r) G d>2A there exists ^ € diA such that -ff(C) < H{i]). 
Let 77 G 5>2.4 and I G the subset of A with smallest cardinality such that a := rj^ E A. The 
definition of d>2A implies |/| > 2; the minimality of / and the construction of S imply r](x) = +1 
and 

ii"(r/^) < H{r]) for all ar G / . (4) 

Now, let X E I, from the minimality of / we get ( := ri^\^^} E ^\A. Hence, = a E A and 
C G \ .4 imply ( E diA. Moreover, from |/| > 2, equation (3) and the property (4) we, finally, 
get H{0 < H{r]). 

Now, by using the same strategy as in Section 3.1 of [6] we get that the minimum of the energy 
on dA is realized in V. Finally, Theorem 1 is a consequence of Propositions 3.4 and 3.7 of [5]. 

In conclusion we can say that in the case of the alternate updating the property (3) implies that 
the mechanism of escape from the metastable phase is not changed with respect to the Glauber 
case; on the other hand one can expect that things could be different if the parallel updating were 
allowed all over the lattice without any restriction [7]. 
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